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Abstract. The transition to turbulence in pipe flow does not follow the scenario familiar from Rayleigh- 
Benard or Taylor-Couette flow since the laminar profile is stable against infinitesimal perturbations for all 
Reynolds numbers. Moreover, even when the flow speed is high enough and the perturbation sufficiently 
strong such that turbulent fiow is established, it can return to the laminar state without any indication of 
the imminent decay. In this parameter range, the lifetimes of perturbations show a sensitive dependence 
on initial conditions and an exponential distribution. The turbulence seems to be supported by three- 
dimensional travelling waves which appear transiently in the flow field. The boundary between laminar 
and turbulent dynamics is formed by the stable manifold of an invariant chaotic state. We will also discuss 
the relation between observations in short, periodically continued domains, and the dynamics in fully 
extended puffs. 



PACS. 47.27.Cn Transition to turbulence 
systems methods 



1 Introduction 



Dynamical system theory and nonlinear dynamics has fig- 
ured prominently in understanding the transition to tur- 
bulence in systems with linear instabilities of the lami- 
nar profile [T] . The various bifurcations and transitions in 
fluids heated from below (Rayleigh-Benard) and in the 
flow between independently rotating cylinders (Taylor- 
Couette) have been studied in considerable details and 
have helped to advance both our understanding of the 
various transition scenarios and of the properties of dy- 
namical systems in general. Pressure driven flow in a pipe 
or the flow between moving parallel plates (plane Couette 
flow) behave differently and are less well understood: the 
linear stability of the laminar profile prevents the imme- 
diate application of the usual bifurcation and transition 
scenarios |2l3l4l5j . So how does the turbulence observed 
in experiments come about? 

Evidence that has emerged in recent years suggests 
that dynamical systems also provide a role model for the 
transition in linearly stable situations [6l5l7j . The applica- 
tion of these ideas to the typically fairly high-dimensional 
situation of pipe flow also raises interesting questions about 
dynamical systems, so that one can again expect a fruitful 
interplay between dynamical systems theory and the fluid 
dynamics of the transition to turbulence. The focus of this 
contribution will be to survey recent progress in our under- 
standing of pipe flow. We will go beyond previous studies 
|8I4I5| by including steps towards a characterization of the 
spatial variability of the system. 



47.27. nf Flows in pipes and nozzles - 47.10.Fg Dynamical 



The main features of turbulence transition in a pipe 
were documented long ago by Osborne Reynolds |9|10j . 
He noticed the possibility to delay the transition by sup- 
pressing perturbations in the inflow region, the intermit- 
tent character of the transition, and the presence of vor- 
tices in the turbulent regions. He noticed that the flow can 
be characterized by a dimensionless combination of mean 
downstream velocity U, radius of the pipe R and kine- 
matic viscosity v, nowadays known as the Reynolds num- 
ber Re = 2UR/v. The intermittent character of the tran- 
sition and the absence of a linear instability of the lami- 
nar profile make the determination of a 'critical Reynolds 
number' above which turbulence can be observed a tricky 
business, as attested by the huge variability of critical 
numbers that can be found in the literature. Reynolds ar- 
rives at about 1800, Prandtl [11] speculates that it should 
be above 1000, MuUin and Darbyshire T7] find long lived 
turbulent trajectories above about 1650, and most refer- 
ence books and textbooks quote numbers in the range of 
2000 to 3000. The operational definition would be that 
the critical Reynolds number is defined such that suffi- 
ciently strong perturbations can induce turbulence if the 
critical value is exceeded. But even with that definition 
there remain problems connected with the transience of 
the turbulent state, as discussed in section 5. 



An interesting feature of turbulence in pipe flow is that 
below Reynolds numbers of about 2700, it does not ex- 
tend throughout the pipe but remains localized in tur- 
bulent puffs of about 60i? length, |13|14] . While some 
progress has been made in describing localized structures 
anf frontal dynamics in various pattern forming systems 
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Fig. 1. Three-dimensional flow inside a puff at Re = 1825 
moving from left to right. The snapshots are separated by time 
intervals of t = AQR/ud and show isosurfaces of the down- 
stream velocity relative to the laminar profile at levels ±Q.lUci 
(negative in blue and positive in red) . The plots are compressed 
by a factor of 5 in the axial direction to show the whole com- 
putational domain of L = lOOi?. Three more frames not 
included for size reasons. 

[15] , the case of localized turbulence is more complicated 
because of the dynamic nature of the localized state [16] . 
Some progress has recently been made for the case of lo- 
calized turbulent stripes in Taylor-Couette and plane Cou- 
ette flow, |17I18| . 

For our presentation below we will draw on some ex- 
perimental observations, and on numerical simulations of 
the full Navier-Stokes equation in domains of length lOOR 
(for the spatially extended puffs) and lOR (for the statisti- 
cal studies in smaller domains), where R is the radius. The 
numerical code has been verified by reproducing the spec- 
trum of the linearized equation, and the turbulent prop- 
erties at somewhat higher Reynolds. The coherent struc- 
tures have also been confirmed by independent studies. 

The outline of the paper is as follows. In section 2 
we discuss the dynamics of a localized turbulent puff and 
present evidence that the interior dynamics can be cap- 
tured by studies in smaller domains. In section 3 we dis- 
cuss coherent states and their bifurcations in domains of 
length lOR. In sccion 4 we present lifetime studies in the 
small domain and relate them to experiment and other 
investigations. In section 5 we discuss the boundary be- 
tween the laminar and turbulent regions in state space 
and the evidence for the presence of an edge of chaos as 
a generalization to the usual basin boundaries. In section 
6 we summarize and give a brief outlook how this relates 
to studies in other systems. 

2 Puffs 

In the Reynolds number range up to about 2700, the tran- 
sition to turbulence in pipe flow shows a striking spatial 
localization [13|14j : the turbulence is concentrated in re- 
gions of about 60R length, which move uniformly through 
the pipe |13|14|1 9 20 21 22 23J. In Fig. [l] we present four 
snapshots of a puff and in Fig. [2] the axial position of its 
center, demonstrating the very uniform speed with which 
it moves downstream. The position of the puff is defined 
via the center of the turbulent in-plane kinetic energy dis- 
tribution. This criterion is very robust as compared to one 
based on a jump in the downstream centerline velocity (see 
[24] ) . and gives almost noise- free position information and 
very accurate translation velocities, as evident from Fig.[2l 
where no smoothing or linear interpolation was applied. 

The axial extension of the structure remains fairly con- 
stant and its traveling velocity decreases form 0.965 U of 
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Fig. 2. Uniform motion of the puff along the pipe. The yel- 
low line is the position of the center of turbulent energy in 
the in-plane velocity components, the dashed line a linear ap- 
proximation to it. Shown are two passages through the pipe 
(which in the numerical simulations is periodically continued 
along the axis). The Reynolds immber is Re = 1800. 

the mean downstream velocity at Re=1800 to 0.94 U at 
Rc=1900. Thus, the turbulent patch is slightly slower than 
the bulk velocity of the fluid. Consequently, when viewed 
from a frame of reference comoving with the puff, there is 
a net flow of fluid from the trailing to the leading edge of 
the puff. 

The velocity field inside the puff can be divided up 
into three regions, a fairly homogeneous turbulent inte- 
rior bracketed by the upstream and downstream front re- 
gions. In order to highlight the different characteristics of 
the velocity fields in these domains, we apply the corre- 
lation function measures introduced in |26j to extract co- 
herent structures. Since in the situation of a localized puff 
the system is not only inhomogeneous in wall-normal but 
also in axial direction we only take out the mean velocity 
as obtained from averages over the azimuthal direction, 
the only homogeneous direction left. Thus, the correla- 
tion functions are evaluated for the local deviations from 
the downstream speed as calculated along a ring at radius 
r = 0.81: 

W2((/), z) = u^(r = 0.81, (f>, z) - {u^{r = 0.81, ^, z))^ (1) 
and 

C{(t),z) = {u,{^j + (j),z)u,{^P,z))^ (2) 

Snapshots of this correlation function in Fig. [3] clearly 
show the three different regions. Intriguingly, the upstream 
and downstream regions are very often dominated by pro- 
nounced coherent structures. On the upstream side, struc- 
tures with three and four streaks dominate. Thus, the cen- 
ter panel of Fig. [3] shows a typical flow configuration. On 
the downstream side, the structures are somewhat shorter 
and occur not quite as often as on the upstream side. 
The dynamics of these coherent structures and their sig- 
nificance for the transition to turbulence was studied ex- 
tensively by Gas van Doorne in his PhD thesis, [27]. The 
central region is spatially and temporally more disordered. 

The disorder in the interior domain is also reflected in 
rapidly decaying spatial auto-correlation functions of the 
three velocity components, see Fig. |4j The downstream 
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Fig. 3. Coherent structures in a turbulent puff at three dif- 
ferent times for Re = 1800. The abscissa runs from z = lOi? 
to 2 = 60i? (half of the computational domain of L = lOOi?) 
and is translated so that the center of the turbulent intensity 
is kept at z = L/S. The ordinate shows the angle 4> from — tt 
to TT with the level indicated by the horizontal line. White 
shading indicates positive values, dark shading negative values 
of the correlation function C((/>, z). Note the rather disordered 
interior and the elongated coherent structures at the upstream 
and downstream edges of the puff. Two more frames not 
included because of size problems. 
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Fig. 4. Autocorrelation functions for the downstream (Czz), 
azimuthal (C^^) and radial (Crr) velocities along the axis eval- 
uated in the comoving frame of reference. The correlations are 
normalized to one for vanishing axial shift Az. They are based 
on the deviation from the time and area averaged mean profiles, 
resolved along the axis. The correlation functions are averaged 
over 3/4 of the computational domain. They include the entire 
puff, and thus also the elongated coherent structures at the 
up- and downstream boundaries. The time average is based on 
about 1000 independent velocity fields. The presented curves 
are for a radial position of r = 0.73, but not much variation 
with r was observed. 

component shows the longest axial correlation length due 
to the presence of streaks. As measured by the level it 
reaches 1/e, it extends over a width of about 2R. The 
correlation functions for the other components drop to 
1 /e of their maximal value over distances of less than one 
pipe radius. The short axial correlation length suggests 
that useful information about the interior dynamics of the 
flow can be obtained by studying relatively short domains. 
Specifically, for a length of lOR the axial correlations in 
the downstream fluid are down to 30 per cent, but the 
computational advantages are enourmous and allow for 
detailed studies of deterministic ^ and statistical anal- 
yses [26] . 

3 Local dynamics in short domains 

The correlation functions of the previous section suggest 
that it should be possible to study the dynamics in short 



domains, thereby eliminiating much of the large-scale spa- 
tial dynamics. Nevertheless, the state space of the system 
remains fairly high-dimensional, and easily reaches 0(10^) 
active dynamical degrees of freedom after the elimination 
of boundary conditions and incompressibility. 

The coexistence of laminar and turbulent dynamics 
for the same parameter values but different initial condi- 
tions in this system suggests that the transition might be 
associated with a subcritical bifurcation (compare |29j). 
For increasing Reynolds number, the usual scenario would 
then imply the appearance of new states in a saddle-node 
bifurcation, followed by the destabilization of the lami- 
nar profile in a collision with the saddle (as realized, for 
instance, in TC flow in the narrow gap limit, [30 J ■ 

In the case of pipe flow new persistent flow structures 
besides the laminar profile have indeed been identified 
[28I31| . They are of the type of full three-dimensional 
travelling waves, see Fig. [5l Three-dimensionality is re- 
quired, since velocity fields with translational symmetry 
in the downstream direction decay and cannot sustain 
any non-trivial flow state. The first structures identified 
were highly symmetric, and contained several pairs of vor- 
tices. Among these, the structure with the lowest critical 
Reynolds number is the threefold set of vortices, with a 
bifurcation near Re = 1250. It is followed by the bifurca- 
tion of the twofold state at Reynolds number Re — 1350, 
and more at higher Reynolds number. Interestingly, the 
friction factors are sometimes higher than the turbulent 
ones, but they tend to approach the laminar values for 
higher Reynolds number. Fig. [6l Less symmetric struc- 
tures may also be identified, and some of them extend to 
lower Reynolds numbers |32j . 

All these states actually belong to continuous families, 
since they tolerate some stretching and compression for 
Reynolds number above their point of bifurcation, |28l7j . 
However, their optimal wavelengths at the bifurcation of 
2.5R and 4.2R for three and two pairs of vortices, respec- 
tively, [28J, can fit well into the short domains studied in 
the next section. 

The bifurcation scenario in pipe flow, however, differs 
from the usual one in that both states that appear in the 
bifurcation are immediately unstable [7j : it is like a saddle 
node bifurcation in an unstable subspace. This has two 
consequences. As a flrst one we note that the absence of 
stable patterns and the hyperbolic dynamics on the sad- 
dle naturally explain the chaotic dynamics. Nevertheless, 
the coherent structures just discussed can show up tran- 
siently in the dynamics J25], where they can be detected 
by apphcation of the indicator ([2]). Studies on time series 
then show that about 20 per cent of the time the flow is 
dominated by the presence of these structures [26]. 

Another consequence of this deviation from the stan- 
dard saddle-node transition scenario is the fact that the 
stable manifold of the saddle is not of co-dimension one 
and therefore does not subdivide state space into the do- 
mains of attraction of the laminar and the turbulent dy- 
namics. Moreover, the saddle does not collide with the 
laminar profile for any finite Reynolds number: the lami- 
nar profile remains stable against infinitesimal perturba- 
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Fig. 5. Coherent states for pipe flow. Shown are velocity pro- 
files averaged in the downstream direction. The deviation from 
the laminar Hagen-Poseuille profile is indicated by colors. Neg- 
ative velocities are shown in blue, positive areas in red. In-plane 
velocity components are indicated by vectors. The regular ar- 
rangement of high- and low-speed streaks reflects different dis- 
crete rotational symmetry classes. 
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Fig. 6. Coherent structures in a friction factor vs. Reynolds 
number plane for pipe fiow. C2 and C3 indicate states with 
two and three vortex pairs, respectively, and the numbers in 
parenthesis indicate different resolutions (see [28]). The dots 
indicate observations on a turbulent flow, and the upper and 
lower dashed lines are extrapolations of the turbulent and lam- 
inar behaviour, respectively. 

tions for all Reynolds numbers [33] . This raises interesting 
questions about the boundary between laminar and tur- 
bulent dynamics to which we turn in section [5l and it 
opens up the possibility that the turbulent dynamics is 
not persistent but transient, a feature we analyze in the 
next section. 



4 Lifetimes 

In the usual subcritical bifurcation scenario the stable 
node forms the basis for an attracting subset in state 
space, where chaos can set in as a result of further bi- 
furcations. However, in the absence of a stable object, i.e. 
if the bifurcation happens within a subspace that is al- 



ready unstable with respect to perpendicular deviations, 
the persistence of the dynamics in the turbulent region 
has to be investigated separately. Indeed, studies by Brosa 
[34] first suggested the intriguing possibility that turbu- 
lence might be transient. As a first step in this direction 
we study lifetimes of a sample of initial conditions in the 
turbulent region. Numerical [35 36 37^ and experimental 
studies [19.20<21I22I36| show that in this region the life- 
times are exponentially distributed. Asymptotically, for 
long times, the probability to find a turbulent state after 
a time t is given by 

P{t) oc e"(*"*°)/^ for t>to, (3) 

with a characteristic decay time r. The exponential distri- 
bution for P(t) implies that the probability density p{t) — 
d\n P(t)/dt to decay at a time t is constant: the flow 
has no memory and decays at unpredictable moments in 
time. This constant decay probability is strong evidence 
for the formation of a chaotic saddle that supports tran- 
sient turbulent dynamics. The characteristic decay time 
T increases rapidly with Reynolds number, but it does 
not seem to diverge at any finite Re. This is confirmed 
in experiments on long pipes and in numerical simula- 
tions as shown in Fig. [7| [36]. Other experimental stud- 
ies on shorter pipes [19I20I21I22] and simulations for five 
Reynolds numbers [38] have been interpreted using a vari- 
ation r(i?e) ~ a/{Rec — Re): however, the extracted pa- 
rameter a and the critical Reynolds number Rec vary over 
a wide range, leaving questions whether the asymptotic 
regime has indeed been reached. Moreover, a reanalysis 
of the data 39J shows that the numerical simulations are 
fairly close to the experimental results from Hof et al [36] . 
The comparison between the data suggests the interest- 
ing possibility that a global bifurcation near Reynolds 
numbers of about 1850 (see also [40]) might result in a 
quantitative change of the scaling function parameters. 
Perhaps previously separated parts of the turbulence sup- 
porting scaffold become connected to raise the lifetimes, 
but without actually turning the turbulent saddle into an 
attractor. 

The observation of transient turbulence complicates 
the identification of a critical Reynolds number for the 
transition to turbulence: ideally, this would be the Reynolds 
number where the flow becomes persistently turbulence. 
The absence of a divergence in the lifetimes rules this def- 
inition out. Alternatively, one might require that more 
than half of all initial conditions remain turbulent up to 
some time: then the critical Reynolds number depends 
on the probability and the observation time. Fortunately, 
these dependences are fairly weak on account of the rapid 
increase with Reynolds number: with the data from [36], 
the characteristic time t (measured in units of 2R/U) 
reaches 1000 for a Reynolds number of 1924, and 2000 
already for 1944. 

5 Edge of chaos 

The coexistence of a chaotic saddle carrying the turbulent 
dynamics and the linearly stable fixed point surounded by 
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Fig. 7. Lifetime distributions for pipe flow, and variation of 
mean lifetime with Re for a pipe segment of length 10/?. Top 
panel: Probability to still be turbulent after some time t for sev- 
eral Reynolds numbers in a semi- logarithmic plot. The straight 
lines indicate exponential fits to the tails of the distributions. 
Bottom panel: Inverse characteristic lifetimes r extracted from 
the exponential fits shown in the top panel. Times in this plot 
are measured in units of 2R/U with R the pipe radius and U 
the mean downstream velocity. 

its basin of attraction naturally suggests to look for the 
boundary between the two objects. In the case of coex- 
isting attractors, this would be the basin boundary. How- 
ever, in order to cover the case of a transient dynamics, the 
concept of the edge of chaos was introduced |41I42I37B5] : 
when crossing from the laminar to the turbulent side by 
increasing the amplitude of the perturbation one notes 
and smooth increase in lifetimes on the laminar side and 
a sensitive dependence on initial conditions on the tur- 
bulent side [35j . Inbetween is a first point where the life- 
time becomes infinite. Because of the chaotic variations 
on the turbulent side, this point sits on the edge of chaos. 
The location of this edge of chaos in state space can be 
expected to be fairly complex: folds have been identified 
[37] and fractal properties cannot be ruled out |44|45|43] . 
Interestingly, different points on the edge of chaos can be 
dynamically connected, and it makes sense to probe for 
the dynamics on the edge of chaos. The evolution of the 
system then shows that there is a relative attractor: it is 
attracting for initial conditions in the edge, but globally 
unstable because of the tendency to either swing up to the 
turbulent dynamics or to return to the laminar one. We 
call this relative attractor the edge state, as it is embedded 
in the edge of chaos. As mentioned before, it extends the 
concept of basin boundary for attractors to the situation 
with an attractor and a saddle, see [33] for further discus- 
sion. The boundary and its dynamics is also relevant for 
determining the smallest perturbation sufficient to trigger 
turbulence, e.g. [46] . 

On the practical side, the boundary can be tracked 
by straddling it with a pair of trajectories where one re- 
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Fig. 8. Visualization of the edge state tracking algorithm for 
pipe fiow. The edge state is bracketed by a pair of initial con- 
ditions in which one initial condition decays and the other 
becomes fully turbulent. After about 200 time units the ap- 
proximation is refined by choosing a new pair of approximating 
trajectories. The temporal modulations in energy indicate that 
the edge state maintains some dynamics and does not become 
a travelling wave. 




Fig. 9. Two snapshots of the edge state at Re — 2160. The 
snapshots clearly show the oscillating motion of the center vor- 
tices between the two high-speed regions near the walls. 

turns to the laminar state directly, and the other swings 
up to the turbulent state. Since the dynamics in the edge 
is unstable, the two trajectories will separate. If the sep- 
aration becomes large enough, a new pair of trajectories 
can be found by refinement, see Fig. [51 This way the pair 
will always stay close to the edge and will approximate a 
trajectory the never leaves it. The velocity fields for these 
trajectories become progressively simpler as time goes on 
|42I37| . In the end, a state consisting of two high speed 
streaks, a low speed streak in the middle, and a pair of 
vortices that drives the streaks, is obtained, see Fig. [9l In- 
terestingly, the two vortices in the middle are dynamically 
active, but do not leave this region. This type of state can 
be followed up to higher Reynolds numbers. 



6 Conclusions 

The studies on the transition to turbulence in pipe fiow 
have confirmed a few expectations (like a subcritical tran- 
sition scenario), but have also revealed unexpected prop- 
erties: the lack of stability of the coherent structures, and 
their transient appearance in turbulent flows, the non- 
persistent nature of the turbulent state, the complicated 
boundary between laminar and turbulent, etc. The ef- 
fort invested into understanding this transition can also 
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be expected to be helpful in other situations: plane Cou- 
ette flow is perhaps the next flow that comes to mind, as 
it shares with pipe flow the linear stability of the lami- 
nar profile to infinite Reynolds numbers. Pressure driven 
plane Poiseuille flow does have a linear instability, but at 
a Reynolds number above the one where turbulence ap- 
pears. All these flows share many similarities, and can be 
profitably studied with the concepts presented here. 

Financial support by the Deutsche Forschungsgemein- 
schaft is gratefully acknowledged. 
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